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It was shown long ago by T. V. Ruzmaikina and A. A. Ruzmaikin in 1 that 
within the framework of a homogeneous and isotropic cosmological model quadratic 
corrections of the gravitational field cannot provide solutions that are both regular 
■ initially and go over to Friedmann type at later times. We find here, by applying 

a dynamical systems approach, 3 the general form of the solution to this class of 
models in the neighborhood of the initial singularity under the above conditions. 
Our starting point in this brief paper is the Lagrangian of the general quadratic 
^v^j | gravity theory given in the form a 

q ■ L(R) =R+ BR 2 + C 7," •'/,', , + /;/,"" '" /,',,,,„,. (1) 

Since for a general spacetime we have the following identity, 
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5 J{R 2 - 1 /,"•'/.'., + R ijkm Ri jkm )^dn = 0, (2) 

in the derivation of the field equations through variation of the action associated 
with (1), only terms up to R^ Ru will matter. If we further restrict ourselves to 
isotropic spacetimes we have a second identity of the form 

■ _ 

5 (R 2 - 3R ij R tJ )V^gdn = , (3) 



which enables us to keep terms only up to R in (1). In the specific model treated 
' herein, we consider a spatially flat universe with metric given by 



ds 2 = dt 2 -b(t) 2 (dx 2 + dy 2 + dz 2 ), (4) 

assumed to be radiation dominated, i.e. P = p/3. Under these assumptions, the 
variation of the action functional constructed using (1) gives the following higher 
order field equations: 
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where k = 6B + 2C. Using (3), (4) the tt-component of (5) can be written in the 
form 
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a the conventions for the metric and the Riemann tensor are those of. 2 
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where 61 is a constant defined from ^|p£ = W (VjT i0 = 0). Note that the 
Fricdmann solution \J2b\t satisfies the above equation. 

Assuming that Eq. (6) has a solution with a regular minimum at t = to, {bo = 
b(t = i ) = and fe = b(t = to) 7^ 0) we can expand this solution as a Taylor series 

K*) = &o + y (t-to) 2 + ^(t-t Q f + --- . (7) 

Direct substitution to Eq. (6) restricts the value of the constant n to n = 
(bi/(bobo)) 2 > b . For this value of k, Ruzmaikina-Ruzmaikin conclude that the 
asymptotic form of the solution to (6) is b(t) w exp((t — t ) 2 /12k) which is obvi- 
ously not approaching the corresponding Fricdmann solution as t tends to infinity. 

We now move on to perform a local dynamical systems analysis in order to find 
the general behaviour of the solutions to Eq. (6) near the initial singularity. This 
analysis is based on the use of the method of asymptotic splittings expounded in 
Ref. 3 As a first step, setting b = x, b = y and b = z, Eq. (6) can be written as a 
dynamical system of the form x = f(x): 

y b\ yz z 2 3y 3 
X = V = ^ Z= 2^-2^-- + 2^ + 2^- (8) 

If a = (a, (3, 7), and p = (p, q, r) we denote by x(t) the solution 

x(r) = ar p = (ar p , (3r q , 7 r r ) (9) 

and by direct substitution to our system (8) we look for the possible scale invariant 
solutions . From all possible combinations, the most interesting is the one with 
dominant part given by 

z 2 zy 3y 3 



and subdominant part 



2nyx 2 ' 2k 

where f = + f sub . The dominant balance (of order 3) turns out to be 



f sub =(0,0,--^ + f I, (11) 



(a,p)=((a,|,-|), (i -\ 
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(12) 



where a is an arbitrary constant. 

The Kowalevskaya exponents for this decomposition, eigenvalues of the ma- 
trix JC = -Df(a) — diag(p), are {—1,0,3/2} with corresponding eigenvectors 
{(4, —2, 3), (4, 2, —1), (1, 2, 2)}. The arbitrariness coming from the coefficient a in 
the dominant balance reflects the fact that one of the dominant exponents is zero 



b Due to a sign mistake, Ruzmaikina-Ruzmaikin in 1 conclude the opposite. 
C A vector field f is called scale invariant if f(ar p ) = r p_1 f(a) 
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with multiplicity one. According to the method of asymptotic splittings, 3 we proceed 
to construct series expansions which are local solutions around movable singulari- 
ties. In our particular problem the expansion around the singularity turns out to 
be a Puiseux series of the form 

oo oo oo 

x(t) =5Z Cli (* - *°)* + *' y ^ =X] c 2i(*-*o)* - *, z{t) = Y,c^(t-t )i-i, (13) 

i=0 i=0 i=0 

where to is arbitrary and c w — a, c 2 o = a/2, c 30 = —a/4. For these series expansions 
to be valid the compatibility condition 

(-2ci 3 + C 2 3 \ 
-c 23 + c 33 = (14) 

-|ci 3 + |c 23 - c 33 / 

must be satisfied. Substitution of Eq. (13) into Eq. (8) leads to recursion relations 
that determine the unknowns cu,c 2 i,C3i. After verifying that Eq. (14) is indeed 
true, we write the final series expansion corresponding to the balance (12). It is: 

x(t) =a(t- to)* + ci 3 (t - t Q ) 2 + a 2 ~ K ^ (t - to)? + ■■■ ■ (15) 

The series expansions for y(t) and z(t) are given by the first and second time 
derivatives of the above expressions. 

Our series (15) has three arbitrary constants and is therefore a local expansion of 
the general solution around the movable singularity to- Also since the leading order 
coefficients can be taken to be real, by a theorem of Goricly-Hydc, 4 we conclude 
that there is an open set of real initial conditions for which the general solution 
blows up at the (finite time) initial singularity at t - Finally, we observe that near 
the initial singularity, the flat, radiation solutions of the higher order gravity theory 
considered here are Friedmann-like regardless of the sign of the R 2 coefficient, while 
away from the singularity they strongly diverge from such forms. 
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